Avalanches Injecting Flux into the Central Hole of a Superconducting MgB 2 Ring 
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Magneto-optical imaging was used to observe dendritic flux avalanches connecting the outer and 
inner edges of a ring-shaped superconducting MgB2 film. Such avalanches create heated channels 
across the entire width of the ring, and inject large amounts of flux into the central hole. By 
measuring the injected flux and the corresponding reduction of current, which is typically 15 %, we 
estimate the maximum temperature in the channel to be 100 K, and the duration of the process to 
be on the order of a microsecond. Flux creep simulations reproduce all the observed features in the 
current density before and after injection events. 

PACS numbers: 74.70.Ad, 74.25.Qt, 74.25.Ha, 74.78. Db, 68.60.Dv 



I. INTRODUCTION 

In recent years it has become clear that instabilities 
are commonplace when magnetic flux penetrates super- 
conductors. In particular, magnetic flux avalanches can 
abruptly penetrate thin superconducting samples in the 
form of tree or fingerlike patterns known as dendrites. 
Such avalanches were first observed in Yl^CusOa;, 
where they must be triggered by a sudden point- like hcat- 
pulse^. Since then spontaneous dendritic avalanches 
have been observed 3 ^^ in Nb, Nb 3 Sn, YNi 2 B 2 C and 
NbN, as well as in patterned Pb thin films^. The re- 
cent interest in the phenomenon was largely triggered 
by the discovery^ that in MgB2 films the avalanches are 
ubiquitous below a threshold temperature T t h ~ 10 K. 
The dendrites disrupt electrical current flow and limit 
the overall current capacity of superconductors^ and 
are thus harmful for prospective applications. On the 
other hand they also represent fertile ground for funda- 
mental research as there are many interesting challenges 
in understanding their nucleation and growth. 

The dendritic instability is now believed to be of 
thcrmo-magnctic origin. Motion of flux in superconduc- 
tors releases heat and increases the temperature. Flux 
motion is facilitated by the increase in temperature, so 
more flux moves which then leads to even more heat- 
ing. If the released heat is not transported away fast 
enough an instability develops whereby large amount of 
magnetic flux rushes into the sample. This mechanism 
was originally invoked to explain global flux jumps where 
the entire sample is filled with a sudden burst of mag- 
netic flux and all magnetisation is los t 11 ' 12 : 13 . However, 
the mechanism has recently also been shown to give lo- 
calised, fingering instabilities^^ Indeed, the instability 
threshold field predicted from this model was recently^ 
shown to agree quantitatively with measurements in both 
Nb and MgB 2 films. 



The dynamics of dendritic avalanches is extremely 
challenging to probe experimentally because of the short 
times involved. Typically the entire process of nucle- 
ation and growth of a single flux dendrite lasts a few 
hundred nano-seconds as measured using a sophisticated 
magneto-optical (MO) set-upii 2 -. It is particularly dif- 
ficult to measure the expected increase in temperature 
associated with the avalanches. On the other hand, it 
is possible to estimate the temperature in the dendrite 
indirectly by designing a suitable geometry where den- 
dritic avalanches can leave a long-lasting, temperature 
dependent imprint. In particular, in superconducting 
rings dendrites may span the entire ring width and thus 
create a heated channel that connects the inner and outer 
edges. During the brief time that this channel is open it 
injects flux into the otherwise screened central hole in a 
process which may be described as magnetic perforation. 
This process depends on the temperature in the channel, 
which may then be inferred from measurements of the 
flux contained in the hole. 

In this paper we report a detailed study on the perfo- 
ration process in a superconducting MgB 2 ring using MO 
imaging. The technique, described in section|TTl provides 
us with a detailed map of the flux distribution over a 
sample- wide area. In section IIIII we present our experi- 
mental results, showing the evolution of the flux $ r<r 
contained within the screened region r < r p during field 
ramps at different temperature, as well as the current dis- 
tribution in the ring. <& r <r p displays abrupt jumps when 
dendrites cross the ring. From a simple model presented 
in section IIVI which connects the temperature evolution 
in the heated channel and the evolution of the current 
during perforations we estimate the peak temperature 
during flux injection. We find that this temperature is 
2.5T C , and that the duration of the flux injection process 
is on the order of one /lis. Finally in section [V] we propose 
a numerical model to describe how the perforation affects 
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the current distribution in the ring. 



II. EXPERIMENTAL SET-UP 

A MgB2 ring was prepared by pulsed laser deposition 
(PLD) on a sapphire substrate. Details of the manufac- 
turing process and typical properties of the superconduc- 
tor can be found elsewher e 17 ' 18 . Using photolithography 
the superconducting film of thickness d = 0.5 /im was 
patterned into a ring with inner radius = 1.10 mm 
and outer radius r\ = 2 mm. 

The magneto-optical indicator was a ferrite garnet film 
(FGF) of thickness 5 /im grown on a gadolinium gallium 
garnet (GGG) substrate and covered with an Al mirror. 
Monodisperse beads (diameter 3.5 /urn) were sprinkled 
on the sample to ensure a small and uniform distance 
between the superconductor and the Al mirror prior to 
positioning of the indicator. This precaution is neces- 
sary since metallic layers are known to suppress dendritic 
avalanche o 18 ' 19 in MgB2. The entire stack was mounted 
in an optical cryostat (a liquid helium flow cryostat from 
Oxford Instruments) and cooled to low temperature for 
observation in a polarisation microscope. We ensured 
sensitivity to field direction by setting the polariser - 
analyser angle a a few degrees away from full extinc- 
tion. Images were captured with a Qlmaging Retiga Exi 
CCD camera with 12-bit pixel depth, controlled from 
a computer. The coil current was controlled from the 
same computer by means of a DAC card used to gener- 
ate a control voltage to the current source. The computer 
also communicated with an external voltmeter (Hewlett 
Packard 34401 A), which was used to monitor the coil 
current. The experiments were then conducted by slowly 
ramping the applied field from to 25 mT while record- 
ing 50 — 200 images at regular intervals. Each time an 
image was recorded, the coil current was read and written 
to a log file before proceeding with the field ramp. 

We calibrated image brightness to field by capturing 
20 — 30 images at low temperature and a point away 
from the sample, while ramping the applied field from 
to 40 mT. The recorded applied field and the image 
brightness were used to fit a second-degree polynomial 
1(B) = a,2B 2 + aiB + ao relating magnetic field to the im- 
age intensity. This polynomial form follows from Malus 
law for light intensity in a polarisation microscope, 

7 = / sin 2 (6l + a)+/ b (1) 

along with the phenomcnological linear relation 9 on B 
for the Faraday rotation at small angles. 

The calibration was performed pixel- by-pixel, which is 
necessary to account for uneven illumination in the field 
of view. The raw images were then processed in two 
stages: firstly, each image was converted to a field map 
using the obtained values of the fitting coefficients {a^}, 
and secondly the flux $ r<rp was computed by summing 



the field at each pixel contained within the relevant re- 
gion. 

Sensor noise leads to observable noise in the field values 
computed for each pixel, but summing over many thou- 
sand pixels to obtain the flux in a region effectively re- 
duces random noise to insignificant levels. A more serious 
concern in quantitative magneto-optical measurements is 
systematic errors (see ref l2Clh . particularly the existence 
of in-plane domains in indicator films which show up as 
triangular features of brighter and darker areas. The do- 
main boundaries move easily in response to changes in 
temperature and applied field, and therefore result in a 
changing background. Wc minimised their influence by 
mounting the indicator very carefully on the sample to 
avoid mechanical stresses, and by applying a weak in- 
plane field. Despite these precautions domain splitting 
did occur in the experiments presented here. This was 
observed as dark regions slowly moving into the field of 
view during the field ramps, thus leading us to calculate 
a too low field value in the affected pixels. Fortunately, 
visual inspection of the images reveals that only a small 
part of the interesting region is affected, never more than 
« 10 % of the total area. 



III. EXPERIMENTAL RESULTS 

Figure [T] (a) shows an example field map recorded at 
T = 6.2 K after applying a perpendicular field to the 
zero-field-cooled sample. Notice the negative field at the 
inner ring edge which is the accumulated stray field from 
the supercurrents. In the lower left quadrant of the ring 
we observe one of the relatively few, but also relatively 
large, dendrites formed at this temperature. In fact, the 
image shown here was the last taken prior to the nuclc- 
ation and growth of a second dendrite which spanned the 
entire ring width and injected large amounts of flux into 
the central hole. This perforation event is illustrated in 
figure [T] (b) showing the difference between the frame in 
(a) and the subsequent frame taken at a slightly larger 
field. The increase in flux inside the central hole, and the 
accompanying reduction near the outer edge, are evident: 
the dendrite brings large amounts of flux from outside the 
sample into the central hole. Remarkably, as dramatic as 
the perforation event is, the flux distribution inside the 
superconducting material is left largely unchanged in the 
process. 

Well below the threshold temperature T t h we observe 
more frequent dendritic avalanches. At T = 3.8 K the 
first dendrite appears at roughly B a = 5 mT. As re- 
ported previously^ the size of individual dendrites in- 
creases with the applied field, until they span the entire 
ring and cause perforations. 

For low applied fields, before the first perforation event 
occurs, the superconductor screens the area contained 
within r pi a characteristic radius where the flux fronts 
from the outer and inner edges meet in a critical state 
superconducting ring. The flux inside this region, $> r <r p , 
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FIG. 1: (left) A field map of the sample at T = 6.2 K and B a = 13 mT. (right) The difference in field between the map shown 
here and the subsequent map in the sequence after a dendrite has perforated the ring and injected a large amount of flux into 
the central hole. A corresponding reduction in the field is seen around the outer edge. In the superconducting material there 
is very little change in the field. The arrows in the field map show the current before and after the perforation calculated from 
the field maps. 



is the sum of the contribution from the supercurrents 
$ se ;/ and the applied field <I>b q = B a A r<rp : 

®r<r p = B a A r<rp + <P self (2) 

Provided no perforations have occurred and the ap- 
plied field is below the field of full penetration B p , this 
sum is exactly zero. Figure [2] shows the evolution of 
$ r<r during ramps of the applied field at four differ- 
ent temperatures. For the highest temperature T = 
7 K (above the threshold T t h for the onset of dendritic 
avalanches in the present experiment) the flux remains 
zero until the applied field reaches B p , after which it 
increases linearly with the applied field. The crossover 
results from the fact that the current in the supercon- 
ductor has saturated and hence the term <& S eif remains 
essentially constant when the applied field is increased 
further. The light gray oblique lines in figure [2] show 
slopes expected from the B a A r<Tp term. Guided by these 
lines one can see that the slope of the measured curve 
is slightly larger than expected if the supercurrents re- 
mained constant. This means there is some relaxation of 
the supcrcurrcnt in the ring as the applied field increases. 

The central observation in this paper is seen in the 
curves obtained at temperatures below Tth- ^r<r is still 
zero up to a certain level, the perforation field, where 
it starts to increase in a step-like fashion, with sudden 
jumps followed by periods of screening. At T = 3.8 K 
there are small and frequent jumps above the perforation 
field, which in this case is w 6.5 mT. The average slope 
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FIG. 2: Measured flux within r p as a function of applied field 
showing crossover from complete screening to a linear increase 
following the applied field. The flux jumps are small and fre- 
quent at low temperature, becoming rarer and larger when 
the temperature approaches Tth- The oblique grid lines have 
slope A r<Tp , with the line through the origin representing the 
contribution to 4 > r <r p from the applied field. The contribu- 
tion to the flux from the ring currents, <&seif, is illustrated for 
one of the curves. The solid black line is the result of a sim- 
ulation with a Kim model for the critical current (see section 
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FIG. 3: (top) The evolution of the current along the radial 
profile and for the same experiment shown in figure [T] (bot- 
tom) The total ring current obtained by integrating the cur- 
rent density across the ring width 



of the flux curve corresponds to the flux increase due 
to the applied field, but $ S ei/ is much smaller than for 
the normal critical state above T t h- This means that the 
avalanches reduce the maximum total current which can 
flow in the ring. 

The curves at T = 6.2 K and 6.8 K show similar be- 
haviour, but with larger and fewer jumps and a slightly 
larger perforation field. In the latter case the jump occurs 
well after the field of full penetration has been reached. 

Clearly the flux jumps are accompanied by an abrupt 
reduction of the total ring current, with a recovery taking 
place during the screening phases. However, it is not only 
the total current which is affected, but also the current 
distribution. Using an efficient inversion scheme based 
on the Biot-Savart law which takes advantage of the fast 
Fourier transform^, we obtain 2-D current distributions 
from our measured magnetic field maps. Superimposed 
on the field map in figure [T] (a) are two sets of arrows 
showing the direction and density of current just prior 
to and just after a perforation. In the former case the 
supercurrent is qualitatively what we expect in a critical 
state ring, with a saturated current in the flux penetrated 
part of the sample and a smaller screening current in 
the Meissner state region. The perforation profoundly 
changes this distribution. Near the outer edge there 
is a significant reduction in the current. Further from 
this edge, but within the flux penetrated region there is 
no discernible change in the current, while it is reduced 
throughout the Meissner region. However, the greatest 
change occurs close to the inner edge, where the current 
actually changes direction. 

All these features are typical after perforation events. 



For the particular image sequence under discussion one 
can see the complete evolution of the current in figure 
[3] Each of the four perforations that occurred in the 
sequence are visible as distinct reductions in the current 
level, with most of the change occurring near the inner 
edge. 

IV. PEAK TEMPERATURE AND DURATION 
OF THE PERFORATIONS 

The experiments give accurate measurements of the 
state just before and just after perforations. Despite the 
fact that the dynamics is far too quick for our experi- 
mental set-up, the fact that we can observe the change in 
&r<r p allows us to make some qualified inferences about 
the details of the process. The starting point is to split 
the process into two phases: (i) the nucleation and growth 
phase when the dendrite tip rapidly moves from the outer 
to the inner edge creating a heated channel, and (ii) the 
flow phase when the resistive channel simultaneously im- 
pedes the supcrcurrents and injects flux into the central 
hole. We emphasise that during the first phase there 
should be little or no change in ^ r <r p . Its duration is only 
determined by the dendrite tip velocity. Ultrafast MO 
measurements^ has been used to measure the tip veloc- 
ity of dendrites in MgB2 yielding a value vu p = 10 5 m/s. 
This tip velocity gives a crossing time 8t\ = 10~ 8 s in 
our sample. 

The second phase is more involved. For a narrow ring 
or a uniform current distribution the evolution of the 
current can be expressed as a balance between the voltage 
drop RI across the heated channel and the induced emf 
voltage LI due to the change in the ring current: 



L- + R(t)I = (3) 

R is an ohmic resistance resulting either from a nor- 
mal state resistivity above the critial temperature T c , or 
a flux flow resistivity originating from moving vortices, 
and is non-zero as long as the channel remains open. The 
inductance L is a property of the ring, and is only de- 
termined by the geometry. The joule heating increases 
the temperature in the channel, a process which can be 
described by the heat equation 

where we have ignored a term describing heat diffu- 
sion into the surrounding superconducting material, c 
is the heat capacity, h is the heat transfer coefficient to 
the substrate, and Tq is the temperature of the substrate 
assumed to be equal to the cold finger temperature. The 
two equations are in fact coupled, since the current in the 
pj 2 - term evolves according to equation Moreover, 
the resistivity is strongly temperature dependent below 
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T c , and hence the evolution of the current is influenced 
by the channel temperature. Physically, one can think of 
the process in the following way: as phase 2 commences 
the pJ 2 term dominates over the heat removal term lead- 
ing to an increase in the temperature. But this has two 
consequences, both of which tend to oppose a further 
increase in temperature: firstly, the heat removal term 
becomes larger for higher temperature, and secondly the 
current is reduced in the ring and hence the heat genera- 
tion is reduced. At some point heat removal equalises the 
heat generation and a maximum is reached, after which 
the temperature decreases until the current density drops 
below J c . At this point the resistivity vanishes, and the 
flux injection halts. 

The fact that c, R and h are all temperature dependent 
means an analytic solution is difficult to find. However, 
equation Q can be used to estimate the maximum tem- 
perature during phase 2 directly by putting dT/dt = 0. 
Following ref \lQ we write the heat transfer coefficient as 
h(T) = ho(T/T c ) 3 and obtain the following equation for 
the maximum temperature T max 



h 



— T max (T max To) 



(5) 



From figure [3] we find that at perforation the aver- 
age current density over the whole ring is J w 0.45 J c = 
2.7 • 10 10 A/m 2 . The heat transfer coefficient was esti- 
mated in referenced^ to be ho — 20 kW/Km 2 on sam- 
ples manufactured in the same way on the same sort of 
substrate as the current sample. If the temperature is 
moderately above T c , the resistivity is independent of 
temperature in this materia l 17 ' 23 i 24 and can be taken as 
p n = IQ- 7 Clm. Finally, with T c = 40 K we obtain 
a peak temperature during the flux injection process of 



105 K. 



We can now estimate the duration of phase 2 from 
equation ([3]). With the above temperature estimate in 
mind R is determined by the normal state resistivity, R = 
p n w/A. It follows that equation ([3]) describes a simple ex- 
ponential decay of current, I(t) = hefore ex p( — pn.w/AL) 
with hefore being the ring current just prior to the per- 
foration. The area A = 45 ■ 10 -11 m 2 is the radial cross 
section of our sample, and w is the width of the resis- 
tive region (the dendrite core). For the jump in flux 
A& se if = {I after ~ hefore)L, we obtain the timescale for 
the flux flow phase as 



st 2 = ^m 1 - 

Wp \1 - self /<&self 



(6) 



The plots in figure [2] can be used to estimate both 
the jump in & se if and the inductance L of the ring. In 
the critical state at B a = B p the ring current density is 
J c everywhere, so the total current is J c dw. Since the 
flux due to the applied field is exactly cancelled by the 
supcrcurrents, we have & S eif = B p A r<r , and we find 
the inductance L = B p A r<rp / J c dw — 4 nH. From the 



first jumps we estimate A$ se z//$ se i/ = 0.2. Notice that 
we now are concerned with the self-flux of the ring and 
consequently the contribution from the applied field must 
be subtracted in the curves in figure [5] The width w of 
the resistive region can be estimated from the core width 
of the dendrites. Barkov et al^- found that typically 
w = 10~ 5 m. With the normal state resistivity used in 
the temperature estimate we obtain 6t 2 « 4 ■ 10" 7 s. We 
stress that this estimate represents a lower bound on the 
true timescale. This is because the temperature is below 
T c during part of the process, with a flux flow resistivity 
Pf which may be orders of magnitude smaller than the 
p n we used to estimate the duration. However, most of 
the current change should occur when the resistivity is 
highest, i.e. p = p n , and our St2 based on a constant 
resistivity is expected to be accurate. 

In ignoring the term in equation ((4]) which describes 
heat diffusion into the surrounding superconducting ma- 
terial, we neglect a possible source of heat loss which 
would lead to a lower maximum temperature. With such 
a term equations © and ^ can be solved numerically 
to yield the complete temperature evolution as well as 
the width of the resistive region with only a single in- 
put from our experiments, namely the change in flux 
<& r<r . However, the calculation requires temperature 
dependent model parameters and detailed initial condi- 
tions for phase 2, which are known only approximately. 
We therefore believe our approximate solution is as reli- 
able as such a full scale model would be. 



V. CURRENT PROFILES DUE TO 
PERFORATIONS 

In the discussion up to this point we have disregarded 
the non-uniform change in the current distribution which 
was pointed out in section [TTT1 In this section we calcu- 
late the current distribution resulting from a perforation 
within the resistive channel model, and find that key fea- 
tures of the observations can be reproduced. The cal- 
culation scheme consists in computing the critical state 
of the ring up to some pre-defined applied field where 
a counter current is injected, reducing the total current 
and changing the current distribution. The trick is sim- 
ilar to the one used in recent work on flux focusing in 
SQUIDs with a narrow sli t 26 i 27 . It is motivated by the 
desire to retain the tractable cylindrical symmetry: the 
counter current provides a simple way to model the effect 
of the resistive slit without sacrificing the symmetry, and 
without considering the dynamical details of the process. 

The counter current distribution should leave the field 
inside the superconducting material unchanged in line 
with the observation from figure[T](b), a constraint which 
is fulfilled by a transport Meissner current for a ring. Of 
course, it is difficult in practice to impose a pure trans- 
port current on a ring: therefore this problem has re- 
ceived no attention in the literature so far, in contrast 
to the analogous problem of a flat strip with a transport 
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currentj 2 ^^. 

The calculations of the field and current distributions 
are inspired by the works of Brand t 31 i 32 ' 33 i 34 , which de- 
tails an efficient numerical method to model the electro- 
dynamics in flat superconductors. The starting point is 
Maxwells equations, written as V 2 A = /10J, and a suit- 
able material law describing the superconducting mate- 
rial, E = E C (J / J c ) n ; E is the electric field, J is the 
current, J c is the critical current of the superconductor, 
E c is a critical electric field, and n is the creep expo- 
nent. While this E — J law describes flux creep, it is 
conveniently used in numerical calculations to describe 
the critical state if the exponent n is high enough. These 
two equations lead to an equation of motion for the cur- 
rent (see refs [32II33I for details) which can be written in 
cylindrical coordinates as 



j{r) = £ Q-HrS) (V)"|^ ~ ^«) dr> (7) 

B a is the applied field. The integral kernel Q(r, r') is 
given by 



Q(r,r') 



r' ' coscj) 



27T J ^/r 2 + r' 2 + 2rr'cost 



--d<t> (8) 



and relates the magnetic vector potential A(r) every- 
where to the current J(r') in the ring via 
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Q(r, r')J{r')dr' 



(9) 



We integrate equation numerically in a straightfor- 
ward manner. The kernel is found by computing equation 
(H]) at discrete points using a non-equidistant grid with 
high grid density near the sample edges^ 3 .. The ensuing 
matrix is then inverted to obtain Q~ 1 (r, r') required in 
equation ([7]) . A technical challenge with this formulation 
is that Q{r,r') diverges for r = r' , but this can be han- 
dled by carefully selecting finite diagonal elements (see 
the appendix for more details). 

The magnetic field distribution is calculated from J(r) 
by first computing the vector potential from equation [5J 
and inserting this into B{r) = dA{r)/dr + A(r)/r + B a . 
The computed flux $ r<rp = J s BdS = 2tt J* B(r')r'dr' 
from a simulation is shown in figure [2] and agrees well 
with the experimental curves. To achieve this close 
matching it is necessary to include a field dependent crit- 
ical current in the simulations. We used a Kim model 
J c = J c o/(l + B/B K ), with B K = 80 mT. This value 
is only a little higher than expected based on previously 
reported 2 ^ magnetisation data on the same sort of sam- 
ples. 

In order to complete the perforation model we also 
need the distribution of a Meissner transport current in 
this cylindrical symmetry. This may be found by invert- 
ing equation @ to find the current J(r) from a magnetic 




FIG. 4: Field (top) and current profiles computed just prior to 
and just after a simulated perforation. The perforation field 
was 0.15 /io Jed, the creep exponent n = 30, and the blocking 
parameter a = 0.13, corresponding to typical values from our 
experiments (see figure [5jl. 



potential A{r). With A(r) = K/r the magnetic field in- 
side the superconducting material is zero, so the desired 
current distribution can be expressed as 



(10) 



It is convenient to choose the constant K such that the 
total transport current equals the total shielding current 
at perforation, i.e. that J Jt r (r)dr = J Ji, e f ore {r)dr. The 
net current just after a perforation can then be written 



Jafter(r) = Jbefore{r) ~ aj tr (r) 



111) 



The parameter a G [0, 1] determines the reduction 
of the total ring current as a result of a perforation. 
This a is analogous to the relative change in flux a' = 
A$ se ;//$ r<rp in the simplified model described in sec- 
tion [TV] Wc artificially inject this current at a predeter- 
mined applied field, and subsequently let the simulation 
proceed normally, i.e. according to equation [7] again. 
Figure 0] shows field and current profiles just before and 
shortly after a perforation. We observe that most of 
the current change occurs near the inner edge. In the 
Meissner region there is some noticeable reduction in the 
current closer to the inner flux front, but towards the 
outer flux front there is hardly any change at all. Near 
the outer edge, though, the current has decreased con- 
siderably. These features compare favourably with the 
experimental current profiles shown in figures [T] and [31 

In figure Owe show computed values of the flux jumps 
size versus the perforation field for different values of the 
blocking parameter a. Also included in the plot is the 
first jump in some experiments. The plot shows that 
the total ring current is reduced by approximately 13 
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FIG. 5: Flux jump size versus perforation field from simula- 
tions and experiments. The filled squares show experimental 
first jumps, and the solid line shows results from simulations 
with a — 0.13. The simulations reproduce the experimental 
fact that the current is reduced by a smaller amount than 



merit elements for the diverging diagonal terms in the 
kernel, equation ([5]). The problem has been dealt with 
in various ways previousl y 31 ' 34 ! 35 ' 36 . One way to obtain 
diagonal elements that work is the following: With an ap- 
propriate grid r, and corresponding weights Wi one can 
integrate Q numerically over the ring width: 



f 1 Q(r,r')dr' wJ^Qyii* 

Jrn 



(Al) 



The right hand side can be split into diagonal and off- 
diagonal terms, and the equation can be rewritten as 



% by perforations. This is considerably smaller than the 
reduction we would infer from the plots of $> r <r p in figure 
[5] using <f> = LI, where the relative change in the flux is 
w 0.2 for the first jumps. This is of course due to the fact 
that most of the current change occurs near the inner 
edge where the effect on the flux in the central hole is 
greatest. The value a = 0.13 also is more in line with 
the plot of total current in figure[3J which shows a relative 
change in current of ~ 0.16 for the first jump. 



VI. CONCLUSION 

We have proposed a method to estimate the tempera- 
ture rise during dendritic avalanches in superconducting 
thin films. This is accomplished using MO imaging to 
measure the flux change in the central hole of a super- 
conducting MgB2 ring caused by avalanches stretching 
from the outer to the inner edges of the ring. Our model 
assumes that the dendrite tip creates a short-lived, resis- 
tive heated channel as it crosses the ring, injecting flux 
into the central hole. From this model we are able to 
estimate the peak temperature during flux injection to 
be T max = 2.5T C = 105 K, the duration of process to be 
5t = 4 ■ 10~ 7 s, and calculate the current profiles before 
and after injection events. 



Acknowledgments 

We thank Paul Leiderer for doing ultrafast MO imag- 
ing on MgB2 samples following our request. AAFO, DS 
and THJ greatfully acknowledge the financial support 
from FUNMAT@UiO and the Research Council of Nor- 
way. 

APPENDIX A: STABILITY OF THE 
NUMERICAL CALCULATIONS 

The stability of the numerical calculations described in 
section [V] hinges on a careful selection of finite replace- 



"')dr' - }2 



(A2) 



The integral of Q(r,r') over r' can be carried out 
by changing the order of integration of <f> and r' . The 
r'-integration can be done analytically; the result is a 
smooth integrand for the (f> integration except for a sin- 
gularity at <j) — 7r, but the integral is well defined over 
the entire interval <f> £ [0, w]. Equation J5]) can be used to 
compute the off-diagonal elements on the right hand side, 
and the resulting set of diagonal elements give stable and 
accurate computations. 

An additional numerical challenge comes from the fact 
that J tr obtained from equation (flQ|) diverges at the 
edges. One way to cope with this for a regular Bean 
model (no field dependent J c ) is to modify the E — J law 
so that for current above |J C | there is an ohmic resistivity 
Pf, in mathematical form that 



E 



PS J 



if J > J c 
'sgn(J) ifJ<J c 



(A3) 



With this small modification the supercritical currents 
quickly relax in the simulations. 

The resistivity does not necessarily have a physical 
meaning in this context since its main purpose is to sta- 
bilise the simulation when the current reaches extremely 
large values. Still, such an ohmic resistivity is realistic 
in superconductors for currents larger than the critical 
current, and the use of the modified E — J law is rea- 
sonable also from a physical point of view. We use a 
Pf/po = 0.001 s _1 . The value will in reality depend on 
magnetic field since it captures the motion of vortices - 
the more vortices moving will generate a larger electric 
field acting against the current. However, for the compu- 
tations the exact value is not crucial provided it is small. 
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